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Abstract 

We study the cosmological evolution of a single BPS D-brane coupled to grav- 
ity in the absence of potential. When such a D-brane moves in the bulk with 
non-vanishing velocity, it tends to slow down to zero velocity via mechanisms like 
gravitational wave leakage to the bulk, losing its kinetic energy to fuel the expansion 
of the universe on the D-brane. If the initial velocity of the D-brane is high enough, 
the universe on the D-brane undergoes a dust-like stage at early times and an ac- 
celeration stage at late times, realising the original Chaplygin gas model. When 
the D-brane velocity is initially zero, the D-branc will always remain fixed at some 
position in the bulk, with the brane tension over the Plank mass squared as a cos- 
mological constant. It is further shown that this kind of fixed brane universe can 
arise as defects from tachyon inflation on a non-BPS D-brane with one dimension 
higher. 

1 Introduction 

In past years, various dark energy models have been proposed to explain the accelerated 
expansion of the universe. Among them, the simplest way is to introduce a non-vanishing 
cosmological constant. However, this model encounters a fine tuning problem in explaining 
the smallness of the cosmological constant. Other than introducing a constant, there are 
dynamical dark energy models, like the Chaplygin gas model, which potentially provides 
a unified description of dark matter and dark energy [U [2] . In this model, the acceleration 
of the universe is explained to be due to the cosmological evolution of some Chaplygin 
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gas whose energy density and the pressure satisfy: pp = —A (A > 0). This energy 
density and pressure relationship can be interpreted as a D-brane moving in the bulk 
in the absence of potential, whose dynamics is described by a Dirac-Born-Infeld (DBI) 
effective action. In the homogeneous and isotropic Friedmann-Robertson- Walker (FRW) 
cosmology, the general continuity equation is d(pa 3 ) = —pd(a 3 ), with a the scale factor. 
From this equation, the evolutionary rule of the energy density of the Chaplygin gas was 
derived: 

where B is an integral constant. Hence, the energy density p oc a -3 at early times 
and p oc \fA at late times. The two stages therefore are suggested to account for dark 
matter and dark energy in cosmology, respectively. However, this model gets negative 
support from observations [U H] because the perturbation spectrum tends to vanish at 
the acceleration stage. 

In this paper, we examine in detail the dynamics of a single BPS D-brane coupled 
to gravity in the absence of potential. We show that the Chaplygin gas can indeed be 
interpreted as such a D-brane slowing down from a high speed in the bulk, with the kinetic 
energy transferred to the internal energy of the brane to fuel the expansion of the universe 
on it. However, in using this interpretation, there is a problem that needs to be solved 
first: what is the mechanism that initially triggers the slow down of the D-brane? In this 
paper, we argue that a mechanism may be the gravitational wave leakage from the brane 
universe to the bulk. 

The investigation on the evolution of a moving D-brane in the absence of potential is 
useful for understanding the dynamics of D-branes that have attained velocities via dif- 
ferent mechanisms in various brane-based cosmological models. In the mirage cosmology 
model j5], the brane attains velocity through interactions with various fields in the bulk 
or on other branes. In the brane inflation model 0, E] , a pair or a stack of D-branes 
and anti-D-branes can roll slowly down the potential, giving rise to inflation, while, in 
the DBI inflation model [91 [10] , the D-branes can roll fast in the bulk. In the tachyon 
inflation model [TTJ [121 EEH HSl US], a single non-BPS D-brane can also attain velocity 
driven by the asymptotic tachyon potential. Our results obtained in this paper indicate 
that the moving D-branes in these models will eventually slow down to zero velocity if 
the potentials become unimportant at the end of inflation. 

Besides the above case for a moving brane, we also find that the D-brane can be kept 
fixed at some position in the bulk if its initial velocity is zero. This case is similar to the 
DGP model [TTJ, [181 [19], m which a three-brane universe is assumed to be embedded in a 
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five-dimensional Minkowski bulk at a fixed position. Here, we verify that a brane indeed 
can remain fixed in the bulk as the universe on it expands. Moreover, this kind of fixed 
D-brane can originate from tachyon inflation on a non-BPS D-brane with one dimension 
higher. 

This paper is organised as follows. In Sec. 2, we present the set-up of the model. In Sec. 
3 and 4, we discuss the two cosmological evolution solutions depending on whether the 
initial velocities of the D-brane is non-zero or zero. The last section gives the conclusion. 



2 The universe on a D-brane 



D3-branes in Type IIB superstring theory are BPS supersymmetry saturated and stable. 
In this paper, we consider the dynamics of a single BPS D3-brane in the 10- dimensional 
bulk, with six extra dimension compactified, in the absence of potential. Concretely, we 
consider the D3-brane to move in a FRW metric in the four uncompactified bulk along 
which the brane extends. Note that the tension of a BPS D3-brane is: T3 = l/[(2n) 3 g s lj], 
where g s is the string coupling and l s = y/a* is the string length. The Plank mass on the 
four uncompactified bulk can be obtained from the 10- dimensional one via dimensional 
reduction [3 [15]: Mp t = V/ (g s l s ) 2 , where V = r®/ (nl®) is the volume of the 6-dimensional 
compactified space and r is the radius of the compactified space. Based on them, we 
introduce the following parameter which is useful for later discussions: 

\2 = T 3 = 9s (2) 

3MJ, 3(27r) 2 /2^' 1 J 



Its magnitude depends only on g s for fixed a' and V. According to [20], a D3-brane affects 
the bulk strongly only within a distance from the brane: r < gl^l s , which is even shorter 
than the string length l s for a small g s or A. In our following discussion, A is indeed 
expected to be small. 

On a D3-brane, we have the low-energy fields including the metric g^ v v = 0,1, 2, 3), 
the U(l) gauge fields and the massless scalars Y 1 (J = 4, 5, • • • ,9) that describe the 
motion of the D3-brane in the six transverse spatial directions. For simplicity, we take 
into account no gauge fields. Denoting the transverse scalars CIS db vector Y = {Y 1 }, we 
can write d^Y ■ d u Y = Y2j d^Y 1 dyY 1 . Since the direction of Y is not important in our 
discussion, we only consider the value d^Yd v Y of d^Y ■ d v Y . Under this consideration, 
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the DBI style action for a D3-brane coupled to gravity is given by 

S = J d^xV^gR ~r 3 J d 4 x^J -det ( 9fiu + d»Yd v Y) 

= TTTT^ / d A x^ \R - lQnGrs^TT^W^y] ■ (3) 

l07TG J l J 

We assume that it is possible to obtain an induced FRW metric on the D3-brane: 
ds 2 = —dt 2 + a 2 (t)dx 2 and as usual Y to depend on time only, being homogeneous and 
isotropic. The Einstein equations lead to the following Friedmann equation 

\a J 6 a z 

The energy density and the pressure on the brane are respectively: 



p= — ^==, p=-rWl-Y^ (5) 
VI -Y 2 

Thus, we have negative pressure on the D-brane for \Y\ < 1 and vanishing pressure for 
|Y| = 1. The pressure and the energy density satisfy pp = — r|, for which D-branes 
moving in the bulk can be viewed as a kind of Chaplygin gas [HE]- The equation of state 
is w = p/p = -(1 - Y 2 ) e [-1, 0]. 

In the flat universe (K = 0), Eq. (jlj) is 

H 2 = -^U, (6) 
yl — Y 2 

where A 2 = t$/ (3Mp t ), as given in Eq. ()2]), and the Plank mass is M Pl = 1/V8^G here. 
We can also have another constraining equation from the equation of motion of Y in the 
FRW metric 

+ 3HY = 0. (7) 



1 -Y 2 

In the next two sections, we will discuss two kinds of solutions to the above two equations 
depending on whether the D-brane is initially moving or fixed in the bulk. 



3 Moving brane: the Chaplygin gas model 

For an expanding universe with positive H, the first solutions derived from Eqs. (j6]) and 
© are 

ln \ tA4 + arctan ( 4 ) = 3A ( t + *o), (8) 



4 



H= ^- r , (9) 

(i -r 2 )i 

where to is a constant. It is clear that H and Y 2 in the solutions are both monotonic 
functions of t, decreasing with the time growing. Specially, when to = 0, Y 2 = 1, H = oo 
at t = and Y 2 = 0, H = A at £ = oo. Therefore, the speed \Y\ of the D3-brane decreases 
while the scale factor a(t) increases as time grows, behaving as the expanding universe. 
The universe on the D3-brane evolves from a pressureless dust-like state with u = to a 
vacuum with u = -1 as time grows. 

Without loss of generality, we set the time to run from t = to t = oo. So the constant 
to decides the initial values of Y(t = 0) and H{t = 0). In the special case of to = 0, the 
initial speed |Y| is equal to the critical value 1 and so the energy density p and the Hubble 
parameter H are both singular at the beginning of time. 

3.1 Dust-like era and late-time acceleration 

Let us take a look at the above solutions at different epochs. At early times, H is large 
H 3> A when t <C 1. In this case, the left side of Eq. ([8]) approximately equals 2X/H and 
the solution can be written 

H~—^ — -. (10) 
3(t + t ) V ; 

Thus, the scale factor and the energy density evolves as, 

2 

a(t) ~ ( 1 + fy « (t + to) 2 ^, p(t) oca~ 3 (11) 

This is the evolutionary rule for dust in the expanding universe. This is understandable 
because the equation of state u ~ as \Y\ ~ 1 at early times. This dust-like phase shares 
a common sense with the pressureless tachyon matter, which forms as the time derivative 
of the tachyon \T\ — >• 1 at the end of tachyon inflation on a non-BPS D3-brane. The 
tachyon matter was investigated as a candidate to cold dark matter [211 1221 E]. 

At late times, the second term on the left-hand side of Eq. (|Sj) is negligible since it 
can only vary within [0, vr/4]. Thus, the solution of H reduces to 

#~ Acoth[3A(t + t )], (12) 

So H decreases to a constant H ~ A at late times, yielding an exponentially expanding 
universe a(t) oc e xt with p oc A 2 . 

Between these two stages, there should be a critical time, at which the universe turns 
from a deceleration stage to an acceleration stage. Let us now determine this critical 
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time. From Eqs. fl6]) and (J7J), we can derive 

a l-^Y 2 

°L = \ 2 2 Y . (13) 
a y/l-Y 2 

It is easy to know that when Y 2 decreases to Y 2 = 2/3, i.e., the time reaches 

0-55 

*<= = — - C ( 14 ) 
the universe enters an acceleration stage from a deceleration stage. The equation of state 
is oj > —1/3 at the deceleration stage (t < t c ) and uj < —1/3 at the acceleration stage 
(t > t c ). 

The above results are consistent with those observed from Eq. (JI| in the original 
Chaplygin gas model [U 12], which is proposed to give a unified description of dark matter 
and dark energy. However, this model has got negative support from observations [3j H]. 

3.2 The trigger of deceleration in the bulk 

In the above analysis as well as in the original Chaplygin gas model [TJ [2], a problem 
is ignored. That is, the reason that the D-brane automatically decelerates in the bulk 
is unclear. If we consider a system including a potential in the bulk, the slow down 
of the whole D-brane may be explained as a result of the brane-bulk or brane-brane 
interactions. In the system concerned here, we may think that a possible mechanism 
could be the gravitational wave leakage from the brane universe to the bulk. As analyzed 
in [23], there indeed are gravitational waves produced in the Chaplygin gas model and 
the evolution of the gravitational waves does not obviously distinguish from some other 
models, like the ACDM model. Since gravity can propagate in the 10-dimensional bulk, 
the gravitational waves produced on the D-brane are possible to flee off the brane, taking 
away part of the kinetic energy and triggering the deceleration of the D-brane. Of course, 
most of the kinetic energy lost due to the deceleration is transferred to fuel the expansion 
of the brane universe such that Eqs. ([6]) and ((7j) are still correct. 

The gravitational waves are associated with the tensor perturbations h^, which are 
symmetric and traceless. Following [2U [22], the gravitational waves evolve as in the 
universe with perfect fluid 

hij + m'h l3 - ^V 2 hij = 0. (15) 

(X 

Defining = a~ 3 / 2 h^exp(ik ■ x)eij, we can rewrite the perturbation equation as 

/ooK=o. (i6) 
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where 

f(Y) = 1°L + *H 2 = ^x/T^. (17) 
JK J 2a 4 4 v ; 

To get propagating gravitational waves, one needs k 2 /a 2 > f(Y). At early times as 
\Y\ ~ 1, there exist gravitational waves for almost all wavelengths. With time growing, 
the scale factor increases and the velocity \Y\ decreases. Hence, the metric perturbations 
with shorter wavelengths tend to freeze out and there are now less gravitational waves 
to leak to the bulk when the velocity of the D-brane decreases. This basically agrees 
with the conjecture given above that gravitational leakage triggers the slow down of the 
D-brane. 

4 Fixed brane universe 

There exist another kind of solutions to Eqs. ([6]) and (J7J). In terms of Eq. ([6]), we can 
write Eq. © as 

Y = -3\Y(1 -Y 2 )f. (18) 

It is clear that, in the special cases \Y\ — and = 1, Y — 0. This means that the speed 
of the D-brane does not change if its velocity is zero or the critical value 1. That is, the 
brane will remain fixed at some position in the bulk if its initial velocity is zero, while it 
will keep moving with the same velocity Y = ±1 if its initial velocity is Y = ±1. In these 
cases, the resulting Hubble parameter H is constant. 

As addressed earlier, the energy density p and the Hubble parameter H on the D-brane 
with non- vanishing tension t% are singular as the D-brane moves at the critical speed ±1. 
Therefore, we only consider the case of Y = 0, in which the action ([3]) reduces to 

S= ltc J rf W^R-6A 2 ). (19) 

It describes a universe with a positive cosmological constant, in the absence of other 
matter fields. This is a special evolutionary mode that was not included in the original 
Chaplygin gas model. 

4.1 Evolution of perturbations 

Although the brane is fixed at some position in the bulk, it can still fluctuate. In what 
follows, we discuss the fluctuations of such a brane universe. We assume that there are 
small perturbations of the form Y = Y(t) + y(t, x) in the homogeneous and isotropic 
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brane universe background. When the brane fluctuates around the fixed position of the 
brane, the action (13]) is 



S 



1 



16nG 



d A X\ 



R-6\ 2 ^l + g^d fl yd u y 



(20) 



The equation of state derived from this action is 



P 

oj = - 
P 



a 2 y 2 



(Vy) : 



(21) 



Correspondingly, there are the scalar metric perturbations in the FRW cosmological 



background: ds 2 = (1 + 2$)dt 2 - (1 - 2$)a 2 (t)df 2 . 
y and $ are respectively 

ij + 3Hy- V 2 y = 
$ + H$ = 0. 



The equations for the perturbations 

0, (22) 

(23) 



Thus, $ oc a 1 drops to zero as the universe expands. Similarly, we define y = a 3 / 2 ^2^y^ exp 
x). The first equation above can be rewritten as 



3 a 3 . 
— + -H' 
2a 4 



y<7 



0. 



(24) 



In the absence of other matter fields in the universe, the scale factor grows as a(t) oc e xt . 
Therefore, this equation reduces to an equation describing propagating waves, whose 
solution for wavelengths q 2 > 9A 2 /4 is 



y(t, x) = a(t) 2 ^^[^4 sin(wt + q- x) + B cos(wt + q ■ x)], 



(25) 



where w = ±y q 2 — 9A 2 /4. Hence, the perturbations y also drops to zero as the universe 
expands. 

From the above solutions of $ and y, we now can determine the evolution of the 
equation of state ( 121 p . At early times of the universe when a is very small, the equation 
of state oj is approximately —2/3. At late times when a is large, the perturbations freeze 
and so the equation of state oj — > — 1. Hence, due to fluctuations of the brane universe, the 
cosmological constant in this brane universe is not always the same everywhere at early 
times. It has a wavy profile of distribution, which leads to oj = —2/3 but not = — 1. The 
equation of state u approaches —1 only at late times. But, in both cases, the expansion 
of the universe accelerates since oj < —1/3. 



S 



4.2 The origin from tachyon inflation 

The Chaplygin gas model and the tachyon condensate model are sometimes combined to 
construct combinatorial model accounting for inflation [261 EH EE]. In this subsection, we 
show that the fixed BPS D3-brane universe discussed above can actually find its origin 
from tachyon inflation on a non-BPS D4-brane. It can be a solitonic object left behind 
the tachyon inflation. 

Let us first take a look at the tachyon condensation process on a non-BPS Dp-brane 
in the Minkowski spacetime [291 ED]- Towards the end of tachyon condensation, some 
solitonic defects form in regions where T = and T = 0. It was verified that these 
defects correspond to some BPS D(p — l)-branes. In other regions, the tachyon field T 
gets infinite as T — > ±1 towards the end of condensation. 

We now show that the situation is similar in the cosmological background. As in 
the usual tachyon inflation model on a non-BPS D3-brane [TTl [T2"l [T3"l [31] . we consider a 
non-BPS D4-brane coupled to gravity 



S T = J d 5 xV=gV(T)y/l + g^d„Td v T, 



(26) 



where V(T) is the tachyon potential which has the maximum value V m at T = and the 
minimum value zero as T — > ±oo. There are two usually adopted potentials derived from 
string theory: V = V m / cosh(/3T) and V = V m exp(—f3 2 T 2 ) with (3 a constant. 

In the FRW metric, the Friedmann equation is H 2 = V/ (3Mp;\/l — T 2 ), with which 
we derive the equation of the tachyon field: 

^T(l — T 2 )3 — (In V)'(l — f 2 ), (27) 



Mi 



pi 



where the prime denotes the derivative with respect to T. For both the tachyon potentials 
V(T) mentioned above, (In V)' = at T = 0. 

Before the inflation starts, the tachyon field fluctuates on the top of the asymptotic 
potential V. As the fluctuation proceeds, we can hopefully find some regions where 
T = 0, i.e., (InV)' = for the above two potentials, and T = 0. In terms of Eq. ( |27j) . 
we know that T = in these regions. Hence, T and T will remain vanishing in these 
regions, where the solitonic BPS D3-branes form. In some sense, the tachyon T is just a 
tachyonic transverse modulus of the non-BPS D4-brane and it will evolve into the massless 
transverse modulus Y of the BPS D3-branes at the end of tachyon condensation. Thus, 
these D3-branes formed at the end of tachyon inflation should be fixed in the bulk with 
Y = T = 0. They provide the origin of the fixed brane universe discussed above. By 
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contrast, in other regions where T and T happen to deviate from zero, T ^ and so T, 
\T\ will grow with time until T 2 — > 1, under which T tends to zero and \T\ tends to cease 
growing. 

5 Conclusions 

We have shown that the transverse scalar Y of a BPS D-brane coupled to gravity is 
stabilised under the condition of the velocity Y = in the FRW metric. The vacuum of 
such a system is in that the velocity and fluctuations of the D-brane are both absolutely 
zero. If the D-brane moves with a non-vanishing velocity in the bulk, it will slow down to 
Y = 0. A mechanism triggering the slow down is the gravitational wave leakage from this 
moving D-brane to the bulk. The kinetic energy lost from the slow down transfers to the 
internal energy of the D-brane, propelling the brane universe to expand. The universe 
undergoes a matter dominated era and an acceleration era as time increases. This is 
consistent with the observation in the original Chplygin gas model. The result indicates 
that the BPS D-brane behaves in an opposite way to a non-BPS D-brane discussed in 
tachyon inflation model [TTl [T2| [T3] : driven by the asymptotic tachyon potential, the 
velocity (the time derivative of the tachyonic scalar) of the non-BPS D-brane in regions 
away from defects increases from to the critical value 1 as time grows. 

On the other hand, if the D-brane has vanishing initial velocity at a position in the 
bulk, it will always be kept fixed at this position. The parameter 6A 2 , which is proportional 
to the string coupling g s , provides a cosmo logical constant. However, in this universe 
with only a cosmological constant, the equation of state evolves from —2/3 to —1 due 
to fluctuations of the D-brane around the fixed location. Interestingly, this kind of fixed 
brane universe is a natural product of defects left behind tachyon inflation on a non-BPS 
D4-brane. 
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